GRAPHENE AND NON-ABELIAN QUANTIZATION 



H. FALOMIR^, J. GAMBOA^'*^, M. LOEWE^-^ AND M. NIETO^ 

Abstract. In this article we employ a simple nonrelativistic model to describe 
the low energy excitation of graphene. The model is based on a deformation 
of the Heisenberg algebra which makes the commutator of momenta propor- 
tional to the pseudo-spin. We solve the Landau problem for the resulting 
Hamiltonian which reduces, in the large mass limit while keeping fixed the 
Fermi velocity, to the usual linear one employed to describe these excitations 
as massless Dirac fermions. This model, extended to negative mass, allows 
to reproduce the leading terms in the low energy expansion of the dispersion 
relation for both nearest and next-to-nearest neighbor interactions. Taking 
into account the contributions of both Dirac points, the resulting Hall conduc- 
■ tivity, evaluated with a (^-function approach, is consistent with the anomalous 

integer quantum Hall effect found in graphene. Moreover, when considered 
in first order perturbation theory, it is shown that the next-to-leading term 
in the interaction between nearest neighbor produces no modifications in the 
spectrum of the model while an electric field perpendicular to the magnetic 
field produces just a rigid shift of this spectrum. 
PACS: 03.65.-W, 81.05.ue, 73.43.-f 



1. Introduction 

Several theoretical ideas have produced a fruitful exchange of interpretations 
and methods between fundamental physics and condensed matter theory, such as 
spontaneously broken symmetry [1] [2l [31 H] or renormalization group methods O 

The recent experimental construction of graphene [10| opens a new connection 
between condensed matter and quantum field theory, since its low energy excitations 
can be represented as massless planar fermions and described by means of a (pseudo) 
relativistic theory. 

Gaphene (See [TTJ [12] and references therein) is a two-dimensional, one-atom- 
thick, allotrope of carbon which has attracted great attention in the last years. 
The carbon atoms are arranged on a honeycomb structure made out of hexagons, a 
structure with a great versatility. Carbon nanotubes, for example, can be obtained 
by rolling the graphene plane along a given direction and reconnecting the carbon 
bonds at the boundaries, giving rise to essentially one-dimensional object. Also, 
the replacement of an hexagon by a pentagon in this lattice introduces a positive 
curvature defect; this allows to wrap-up graphene to give fullerenes, molecules where 
carbon atoms are arranged on spherical structures. Notice that the accumulation of 
graphene layers, weakly coupled by van der Waals forces, constitute the well-known 
graphite. 

This peculiar material has been theoretically predicted by Semenoff in 1984 [13] 
(See also [14j) and experimentally produced in the lab in 2004 jlO] . 
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The electronic properties of graphene are the result of the sp^ hybridization 
between one s and two p orbitals, which leads to a trigonal planar lattice with 
a a bond between carbon atoms separated by 1.42 Angstrom. This filled band 
gives the lattice its robustness. The third p orbital of the carbon atom, oriented 
perpendicularly to the plane of the graphene, gives rise to the tt band through 
covalent bounds with neighboring atoms. Since this p orbital contributes with only 
one electron, the tt band is half-filled in neutral graphene. 

The low energy excitations of graphene accept a description as states of chi- 
ral massless Dirac fermions with a pseudo-relativistic linear dispersion relation, in 
which the speed of light c is replaced by the Fermi velocity, vp ~ 10~^c. Then, the 
Lagrangian describing these low energy states in the presence of an Electromagnetic 
field is similar to that of QED for massless fermions, shearing therefore some of its 
peculiarities. In particular, when a magnetic field is applied perpendicularly to the 
plane of graphene, an anomalous integer quantum Hall effect |15[ 116) takes place, 
which has been experimentally measured [17[ 118) . 

It is not the goal of the present paper to discuss the physics of graphene from 
first principles. Rather, we will consider a simple effective non-relativistic Hamil- 
tonian, suggested by a particular deformation of the Heisenberg algebra (non- 
commutativity of momenta, consistent with the introduction of an external con- 
stant non-Abelian magnetic field), which could be useful to describe the low energy 
excitations produced by the dominant nearest and next-to-nearest neighbor inter- 
actions in graphene. 

The non-commutativity of spacetime is an old idea [19) . the first example of 
which was probably discussed by Landau in 1930 [20]. It has been revived in recent 
years within the context of string theory [21) and since then, non-commutative 
field theories have attracted much attention in various fields such as Mathematics, 
Theoretical Physics , [H |23l|2l|25] , and Phenomenology [26] . 

The breaking of commutativity of the position operators and the representations 
of the algebra of the non-commutative space-time coordinates has been studied in 
[2 7) . The non-commutativity in the momenta algebra we are interested in can be 
related to the deformation quantization of Poissonian structures developed in [5S] 
and considered as a kind of magnetic quantization [29[ 130] . 

All these researches have stimulated the construction of new models in quantum 
mechanics [31) . which have opened new routes to explore, for example in supercon- 
ductivity [6]. Also, a massless Dirac-like Hamiltonian in a generalized noncomniu- 
tative space and its relation with graphene has been considered in [5^ . 

Recently, some models based on a kind of nonstandard deformation of the Heisen- 
berg algebra, which can be realized by shifting the dynamical variables with the 
spin, have been studied in [33l [34]. In the following, we will consider a similar 
deformation, but concentrated in the commutators among momenta, which can be 
interpreted as the introduction of a constant non-Abelian magnetic field [35 . 

In the next Section we present the model and derive the Hamiltonian. In Section 
|3]we study the free case and in Section |4] we introduce a constant magnetic field 
and solve the corresponding Landau problem. In Section [5] we apply our results to 
describe the low energy states associated to the leading nearest and next-to-nearest 
interactions in graphene. 
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In Section [6] we evaluate the associated Hall conductivity employing a C- function 
approach, finding that the result is consistent with the anomalous integer quantum 
Hall effect present in this material. In Section [7] we are also able to show that 
the next-to-leading (quadratic) term in the nearest neighbor interaction does not 
change the spectrum at first order in perturbation theory. In Section |8] we comment 
on the case where crossed electric and magnetic fields are present and, finally, in 
Section |9] we establish our conclusions. 

At the end of the paper. Appendix |X] is dedicated to the Lagrangian of the 
model and studies its symmetry, conserved current and the generating functional 
of Green's functions. It also comments on the weak field and gradient expansion of 
the generating functional, its relation with the Hall conductivity and the topological 
considerations involved. 

2. Deformation of the Heisenberg algebra 

We consider particles leaving on a plane whose dynamical variables satisfy the 
deformed Heisenberg algebra given by 

(2.1) [X,,X.j\=Q, [X„P,]=z%, [P„P,] =2z0%-3a3, i,i = l,2, 

where the momenta commutator is proportional to the pseudospin and is a 
parameter with dimensions of momentum (For convenience, we take h — 1 and 
return to full units when necesary). 

Notice that these particles are described by wave functions with two components, 

= ^ ^ ^ e L2 (M^) (giC^, and these operators have the structure of 2 x 2 matrices 

on C2. 

The deformed algebra in Eq. (|2.1I) can be realized by defining 

(2.2) X,■.= x^(g>l2 Pj := «) I2 + 6* «) cr« , 

with , i = 1, 2 the usual (commutative) coordinates on the plane and pi — —idi, 
operators on L2 (K^), and Ci , i = 1,2 the first two Pauli matrices. We will also 
write cr := (Ils ® ci, ® (T2). For notational convenience, from now on we will 
avoid the explicit indication of the symbol ®, which can lead to no confusion. 

Our aim is to consider the direct generalization of the Hamiltonian of a (non- 
relativistic) particle of charge e and mass to, minimally coupled to an external 
(2+1-dimensional) electromagnetic field, {Aq, A := (Ai, ^2)}, which is constructed 
by means of the replacements Xi ^ Xi, pi ^ Pi. 

For the time being, we make Aq — 0. Then, 

2to 2m 
where the electromagnetic field is taken in the Coulomb gauge, V • A = 0. 
The Hamiltonian can also be written as 

(2.4) (P~eA)^ +VFa-ip-eA) , 

2to 

where we have defined the Fermi velocity 

, ^ Q 

2.5 VF := — 

TO 

and subtracted the constant 9'^ /m. 



4 



H. FALOMIR-*, J. GAMBOA^'^, M. LOEWE^'° AND M. NIETO^ 



In the m — >■ oo limit, with fixed vf, the resulting linear Hamiltonian is appro- 
priate to describe the conducting effective particles in graphene around the Fermi 
points [131 [TTJ [121 US] , which justify our proposal. Notice that this limit does not 
correspond to a small but rather to a large deformation of the commutator [Pi,P2]. 

Notice also that the modification of the Hamiltonian in Eq. (12. 4[) can also be 
interpreted as the introduction of an SU{2) non-Abelian constant and uniform 
"magnetic" field ~ 9^. Indeed, the SU{2) transformations relate both components 
of the wave functions, and the commutator of covariant derivatives gives 

(2.6) [pi - eAi + 0ai,p2 - e^2 + Oa^] = le (9x^2 - ^2^1) + '2.i6'^<Jz ■ 

In this sense, if we take a constant magnetic field B ~ {diA2 — 92^i), the system 
we are considering is a kind of non-Abelian version of the Landau problem [35] . 

In Appendix [A] we describe the Lagrangian of this model, study its symmetry 
and conserved current and discuss the relation of the asymptotic expansion of its 
generating functional with the Hall conductivity. 



3. The free case 

In this Section we consider the free case, with A = 0. The Hamiltonian reduces 



to 

(3.1) 



H = h o- • p . 

2m 



We propose solutions of the form 

(3.2) Vk(x) = e^'^-^xi^) , 
with x(k) G C^, which replaced in 

(3.3) [iJ-£(k)]^k(x) =0 
lead to 

k2 



(3.4) 



2m 



i;F<T.k-£(k) x(k)-0 



Nontrivial solutions require 



(3.5) 

which gives 
(3.6) 



det 



2m 



WFcr •k-5(k) ^ = 0, 



2 1 2 

Vf k 



^ ' 2m 



Then, we get the following dispersion relation (approximately linear for small |k|, 
see Figure 1) 



(3.7) 



£{k) = ^±VF |k| 



which, replaced in Eq. p.4p for k ^ 0, shows that the pseudo-spinor x(k) has 
definite pseudo-helicity. 



(3.8) 



(,T.k)x±(k)=±x±(k). 
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Figure 1. The dispersion relations for the two branches of solu- 
tions for the free case. 

On the other hand, for k = the two linearly independent vectors [ ^ ) and 



are just constants solutions with vanishing eigenvalue. 



The Hamiltonian in Eq. p.ip commutes with the effective angular momentum, 
the generator of a U{1) symmetry, 

(3.9) J -id^ + i fT3 , [H,J]=0. 

Under a rotation on the plane, the wavefunction in Eq. p. 21) changes into 

(3.10) U{^)iP{^) e'^^'^i; (i?(i?)"'x) = e'^'i"^^^^'''') e' ^ "^x(k) . 
For ?9 = 27r we have i?(27r) = I3 and e''^'^-^ ^ go, we get 

(3.11) U{2Tr)TP{x) = -7/>(x) . 
Therefore, these particles can be considered as fermions [TT] . 

4. Constant magnetic field perpendicular to the plane 

In this Section we consider the electromagnetic vector potential of a constant 
magnetic field orthogonal to the plane, 

(4.1) A = Bxie2 ^ 91^2-92^1=5 and V-A = 0. 
In this case the Hamiltonian can be written as 

(4.2) 2mH ~ pi^ + (p2 — eBxif' + 2mvpaiPi + 2mvp(j2 {p2 — eBxi) , 

which clearly commutes with p2 ■ This allows to look for generalized eigenfunctions 
of the form 

(4.3) V(x) = ^ $(xi) , 



where $ = 
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The eigenvalue equation for the Hamiltonian, [H — £)tp — 0, reduces to the pair 
of coupled differential equations 
(4.4) 



pi' + {eBf 



Pi 



-1 




-1 





k 

eB 
k 



with A — 2m£. 

At this point, it is convenient to change the remaining variable xi in favor of 
q y^\eB\{xi 
(|4.4p writes as 



k_ 

eBJ 



So, pi = ^y\eB\p, with p 



, and the system in Eq. 



(4.5) 



{\eB\ [p2 + g2] _ X} ip{q) = -2mvFV\eB\{p + zsgn(eB) q}x{q) , 



{\eB\ [p^ + q^] -X}xiq)^-2rnvFV\^\{P~^sgii{eB)q}ip{q). 



For simplicity, in the following we will take eB > 0. The case eB < can 
be obtained from the previous one by simply interchanging the components of the 
solutions, Lp{q) o x('z): can be easily seen from Eq. ()4.5p . So, we consider 



(4.6) {eB [p' + q^] - \} 



The presence of twice the Hamiltonian of a harmonic oscillator of frequency 
1 in the brackets on the left hand side of Eqs. ()4.6p . operator with eigenfunctions 
4'niq) = e"* ^"^Hniq) {Hn the Hermite polynomials ) and eigenvalues 2ri + l, for n — 
0, 1, . . . , suggests that if{q) ~ 4>n{q) and x{q) ^ 4>n'(,q), for some n, n' . Moreover, 
since the Hermite functions satisfy 

{(kn'iq) + #n(g) = 2r7,(/)„_i((7) 
K'iq) ~ qKiq) = -4>n+i{q) 

one concludes that these solutions are of the form 

(4.8) Lp{q) = ci(j)n+i{y) , X{q) = C2(pn{q) , 

for n = 0, 1, 2, • • • , with Ci, C2 G C. Replaced in Eq. (|4.5p . wc get the homogeneous 
system of algebraic equations 



(4.9) 
with 

(4.10) 



M 



Cl 
C2 



M = 



/ e-B(2n + 3)-A i2mvF\/eB \ 

-i2'mvFVeB{2n + 2) eS(2n+l)-A 



V 



The eigenvalues are determined by the condition 



(4.11) 



detAf = [2eB{n + 1) - A] - [eBf - 8m^VF^eB{n + 1) = , 
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which implies that 



(4.12) 



f — 



2m 



1 r 

z 



with s = ±1 and 
(4.13) 



mvp 



Appropriately normalized, the generalized eigenfunctions are written as 



(4.14) Vfc,n,.(x) 
where 



Kn.s e 2 



V 2z 



Hn+i{q) 



\ 



2-2- 



(4.15) 

and 
(4.16) 



^^(f + v/rT8z2(„ + i))'- 1 



Indeed, taking into account the orthogonality relations for the Hermite functions 
it can be easily verified that 



ilpk,n,s,i^k',n\s') = Sn,n'Ss,s'S{k - k') . 



(4.17) 

Finally, notice that there is another solution $o of Eq. (|4.6p whose components 
are given by 

1/4 



(4.18) 



6o(g), x{q)=0- 



Indeed, since Ho{q) = 1, we have 

(4.19) {p - iq)Mq) = -<dg + q)e-«'/2 - 

and Eq. (14. 5p reduces to 



(4.20) 



eB 1 
'm 2 



+ q'\ - Mq) = 



eB 1 
~m 2 



-£o}Mq) = Q- 



This implies that £q = eB /2m, which is independent of vp 
One can also verify that 

(eBV'^ ( Mq) 



(4.21) 

satisfy 
(4.22) 



V'fe,o(x) 



'2tt V I" 



(V'fe.n.s, V'fe'^o) = , (V'fe.o, V'fc'.o) = (S(fc - A:') , 



As mentioned in Section [51 these results can be interpreted as the solution of a 
non-Abelian version of the Landau problem, in which we have also introduced a 
constant and uniform non-Abelian magnetic field. Moreover, in the m ^ oo limit, 
these eigenfunctions are similar to the solutions found 1 361 for the Dirac equation 
in a constant magnetic background. 
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Notice that one can take appropriate linear combinations of the generahzed eigen- 
functions in Eq. (|4.14p and ()4.2ip so as to construct a manifestly complete set of 
generalized vectors in our space, 



(4.23) 



■4>n{q) 







;fce R,n = 0,1,2,, 



27r ' V 1 

where (pnil) are the Hcrmitc functions. Then, the set of generalized eigenvectors of 
H we found is also complete, which ensures that in our analysis we got the whole 
spectrum of the Hamiltonian. 

In the following, we consider the large-z (large-m) limit we are interested ir0- 

For z ^ 1 the energies are given by 

(n + 1) 



(4.27) £n,s = VF 

and the eigenfunctions reduce to 



eB ^^s^2(^) + + O (z-2) I 

H„+i{q) + O (z-^) 



(4.28) Vfc,n,.(x) 



/2^ ^y/{n + iy. 



2 e 

--v'8(n+l)i/„(g)+0(z-i) 



Negative mass. Since we are interested in the description of low energy states 
and the construction of the solutions of Eq. (14. 4p is independent of the sign of m, 
we can explore the behavior of the system for m < 0. Indeed, it is safe to change 
the sign of m in the solutions, which corresponds to the replacement z — ^ —w, with 
w — \m\vF/VeB > 0, in the last line in Eq. (|4.12|) maintaining vp > 0. We get 
(See Figure 2) 

- (ypVeS) — \n + 1 + -^1 + 8w'^(n + 1) = 
\ / w I 2 



(4.29) 



(vfV^) |-,sV2(n+l) - ^ + O {w-') I 



-'^For z <^ 1 wc have equally spaced levels in each branch, with slightly different slopes, 
(4.24) 



e„,s = —\(n + l + -) +2sz'^in + l) + 0(z'^)] 
For the eigenfunctions wc get 



V'fc,n, + l(x) 



2~2~2 



: e 2 X 



(4.25) 



and 



(4.26) 



/2¥ ^^(" + 1)! 
z^{n + l) + 0{z^)]H„+i(q) 
-2iz{n + 1) + &iz^(n + lf +0 {z*)] H„{q) 



V'fc,n,-i(x) = —e-^x 



I 



\ 



2-K ^Vn\ 
z-3z^in + l) + 0{z'')] H„+i{q) 
11 



1 _ ^2(„ + 1) + — z^{n + 1)2 + O (z^) 



Hnig) 

The 2 — > 0-limit is consistent with the usual Landau levels problem. 
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Figure 2. The energy levels £n.s for negative mass, in units of 
vpVeB and for w — 100. 



On the other hancO, So = -eB/2|TO| < 0. 

5. The relation with graphene 

The structure of graphene [TT] can be seen as a triangular lattice with a basis of 
two atoms per unit cell. The lattice vectors are 

1 \ 3/1 

-1 



(5.1) ai=2'^(j-)' ^2 = ^a _ 

where a « 1.42A is the lattice constant (distance between nearest neighbor carbon 
atoms). The vectors characterizing the reciprocal lattice (whose elementary cell is 
the Brillouin zone) are given by 

The superposed triangular lattices form an hexagonal honeycomb array of carbon 
atoms where the nearest-neighbor vectors are 

The tight-binding model for graphene, where it is assumed that electrons can 
only hop to both nearest ((i,j)) and next-to-nearest (((?,j))) neighbor atoms, is 
described by the Hamiltonian [11] 



(5.4) 



^Notice that £"0 — > for m — > 00. This is consistent with the well known existence of zero 
modes in the spectrum of massless Dirac fermions coupled to gauge fields, discussed in different 
contexts in |37| . |38| and | 39| . for example. 
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where t and t' are the hoping energies, a| ^ and ai^s are the creation and annihilation 
operators of electrons in the site i with spin s belonging to the triangular sublattice 
A, and similarly for the sublattice B. By Fourier transforming these operators 
and diagonalizing the resulting expression, one gets the band structure (dispersion 
relations) [IT] 



(5.5) E±{k)^±t^m-t'[fik)-3] , 
with 

(5.6) /(k) 3 + 4 cos (^^Y^^ cos ^^^20^ + 2 cos (\/3fc2a) > . 
The minima of /(k) in the Brillouin zone correspond to the Dirac points 

where /(K) = = /(K'). Then, the bands corresponding to each sign in the 
dispersion relations, Ec. (|5.5p . touch each other at the Dirac points. Notice that 
these two independent Dirac points have been chosen so that they are related by 
the reflection of k about the /ci-axis (fc2 — ^2) 

A series expansion of E± (k) around k = K leads to 

-a|k| - - a^k^ sin(36l) 



(5.8) EsiK + k) = st 



+ t' 


--a^k^ + 3 




4 



o(ikn. 



^2 ' ' 8 

where tan(^) — k2/ki and s — ±1. 

The expansion around the second Fermi point, K', leads to the same expression 
with k reflected about the fci-axis. This corresponds to the change 9 — > —9, which 
changes the sign of the the second term in the first bracket in the right hand side 
of Eq. jET 



We now turn to the comparison with the model developed in the previous Sec- 
tions. From Eqs. (|5.8p and p.7p . it is seen that we can identify vp = | ai. Moreover, 
the next-to-nearest neighbor contribution can be represented in our (free) model by 
the mass term through the identification — | t'a^ = up to a rigid displacement 
of the spectrum in an energy 3i' (which can be subtracted to restore the particle- 
hole symmetry at linear order in |k|). This means that we must consider a negative 
mass m our model, m = -2/ (gi'a^) < (0. 

On the other hand, the second (quadratic) term in the nearest neighbor contri- 
bution depends on the direction of k. We can not reproduce this behavior within 
the framework of our model since it is rotationally invariant. Rather, we must 



■^According to the values of these parameters reported in [111 , we have 

(5.9) a=1.42A, t = 2.8eV, t' = O.leV . 
Then, 

(5.10) J? = 600 Km/sec, mc^ = -4.3 X 10*^ eV , 
or, in natural units (vp/c — > vp, mc/h m), 

(5.11) vp = 2xW^^, m = -8.53 X 10^^ meter-l . 

This is a rather large value for the (negative) mass of the quasiparticles, since |r7i| ~ 32.9me. For 
example, for an external magnetic field B = 10 Tesla, we get w = ^"^'^ \/ = 1385. 
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treat it as a perturbation on the solutions we found for our model. This will be 
considered later. 

Then, we will assume that the effective Hamiltonian for the low energy states of 
this system around K in the presence of an electromagnetic field, Hj^, is obtained 
though minimal coupling and coincides with the one given in Eq. (j2.4p . 

Since the application of our model to graphene requires w ~ 10'^ (See footnote 
[3]), the large w expansion in Eq. (|4.29p gives a good description of the self-energies. 
So, we get for quasi-particle states 



(5.12) 

and for hole states 



1 72(^^-^+0(^-2) I , 



(5.13) wfV^M -V2(^i+1) + 0{w 



Moreover, there is an additional (hole) state at energy Sq — — "^2^{f^ ^ slightly 
below the zero energy level (see Eq. (14. 20^ *1. Indeed, the difference in energy for 
the two first hole states is v^(A/2 + 0(w-i))g. 

5.1. The second Fermi point contribution. As previously mentioned, the dis- 
persion relations near the Fermi points are related by the reflection of momenta 
about the /ci-axis. Then, for the linear terms in the effective free Hamiltonian 
describing the states around the second Fermi point we can write [llj 

(5.15) Vp {(TlPl - (T2P2) = C* ■ P = -0'2 (vp (T • p) 0-2 • 

Accordingly, we will assume that the effective Hamiltonian describing the quasi- 
particles near the second Fermi point in our model, Hj^i , can be obtained from 
through the transformation 

(5.16) Hk' = -<T2Hk(72 - _1p_£^ + yp (T* (p eA) . 

Zm 

Notice that both iJx and Hk' are left invariants under time-reversal times parity 
transformations, TV (See the discussion in Ref. (12) about the symmetries of this 
system). Indeed [12], 

iTV)HK{TV)^ =<JiHk*(Ji = 

(5.17) ( ( - 1 

= ^1 <^ ^-^^ +VFa*-{p- eA) Ui = Hk 

and similarly 

{TV) Hk' {TV)^ = -a2* [TV) Hk {TV)^ 02* = 

(5.18) 

= — cr2^^K0'2 — H^i . 



^For ui = 10'' and _B = 10 Tcsla, we have in full units 

, vp j eB 

(5.14) vpVeB^hc—J 

c y h 

This leads to So = -2.43 X 10"^ eV and £0,+ = -6.89 X lO'^ eV. 



12 



H. FALOMIR-*, J. GAMBOA^'^, M. LOEWE^'° AND M. NIETO^ 



Evidently, the spectrum of Hk.> is obtained from that of Hk by a reflection about 
the origin and its (generalized) eigenfunction are just (T2 times those of i?K- This 
interchanges quasi-particles with holes. Indeed, 

(5.19) Hk' (cr2V'fc,«,s) = -0-2 (ffK^fe,«,s) = ~'£n,s {(^2'lpk,n,s) ■ 

Moreover, 

(5.20) Hk' (aaV^o) = -<J2 (i^K^o) = So ('^aV'o) , 

with an energy £'q — —Sq ~ ^''^i^ > 0. This state corresponds to a quasi-particle 
of energy slightly above zero. 

Therefore, taking into account the eigenstates of both H-^ and iJx' we get, for 
quasi-particles and for holes, an almost doubly degenerate spectrum, except for one 
state of quasi-particle and one state of hole with energies near zero. For positive 
energy states we have 



(5.21) 



£•„,_ = (vfV^) I V2(n+1) - ^ +0{w-^)^, 
-£n,+ = («F Vis) I v/2(n+l) + + O (w-^) I 



where the gap between contiguous states is 

(5.22) A£n = ^(vFV^) \^^^^ + 0\ 



As we discuss in the next Section, this spectrum reproduces qualitatively the 
anomalous integer quantum Hall effect found in graphene [HI [TBI [H EH] : which 
shows a nonvanishing Hall conductivity for a small (positive or negative) Fermi 
level. 

6. The Hall conductivity 

As discussed in Appendix \^ the Hall conductivity has a topological character 
and can be calculated from the weak field and gradient expansion of the effective 
action of the system. However, since we know the exact eigenvalues of energy of 
our model, we will employ a more direct evaluation method, based on the relation 
established between the conserved current and the external electromagnetic field, 
Ecs. (|XT3l) and ((XTl) . 

For our purposes, it will be sufficient to consider the mean value of the den- 
sity jo{^)- Then, as in the previous Sections, we will take E = and B 
perpendicular to the plane of the system. 

Let us start by writing down the partition function of the particles around one 
Fermi point at inverse temperature (3 and chemical potential fi. It can be obtained 
through a Wick rotation of the generating functional Z[A] (see Eq. (jA.Sp ) to the 
2-1-1 Euclidean space by means of the replacements x''^ = t —i' —it, Aq lA^, 
maintaining unaltered the other coordinates and components of the gauge field. 

Since, in our case, Aq = = Ai, A2 — Bxi, from Eqs. (|A.1|) and (IA.8I) . this 
Wick rotation leads to the partition function 



(6.1) Z{l3,fi,B)-~ Vil^^VtPeJo J I Ct 
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where H is iJK (or ifK') and the functional integral is performed on the set of 
configurations of the ferniionic field which satisfy anti-periodic boundary conditions 
on [0,/3]. Then, 

^^(/3,M,i3)= f dtj£x{rpH) = 

(6.2) ^ " 

= P j d^x-^ MP, ii,B)=p j d^x^B, 

where use has been made of Eq. (|A.14|I . 

Our goal is now to evaluate the partition function as the functional determinant 

(6.3) Z{(i,^jL,B) ^ Dei {D) 

where D ~ — iJ is a differential operator defined on a domain of anti-periodic 

functions of t g [0,/3]. Even though D is not symmetric, since H does not depend 
on T this operator has a complete set of orthogonal generalized eigenfunctions 
constructed as 

(6.4) = e-'^i^ , *i,fc,o = ^ e-*'^'^ Vfc,o , 

with ^ e Z, fc e M, n = 0, 1, 2, . . . , s = ±1, and where the Matsubara frequencies 

(6.5) ~j3 \ ^ 2) ' ^ ^ -^-i-i ■ 
Then, the eigenvalues of D are given by 

(6.6) 

where £n,s and f o are the eigenvalues of H, studied in the previous Sections. 

Notice that Xi^n,s and Aj^o are independent of k. Therefore, in evaluating T)et(D) 
we will forget about the index k and, at the end, take into account the degeneracy it 
introduces per unit area, given by the number of flux quanta per unit area [40| 141) , 
A = eB/2n (or A = eB/h in fuU units). 

We define Det{D) as the (gauge invariant) C-function determinant [HISS], 

(6.7) log Det (£>):= --^Tr- 

du 

where A is an arbitrary mass scale, the trace is evaluated for 5ft(u) sufficiently large 
and u — > stands for the analytic continuation of its derivative to a neighborhood 
of u = 0. In our case, 

Let us first consider the contribution of the second term on the right hand side 
and evaluate 

(— =l-\ A A 




/=-oo ^ ' ;=o ^ ' 1=0 
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where use has been made of the properties of the Matsubara frequencies stated in 
Eq. (|6.5p . We can also write 



(6.10) 



OO 

E 

^— — oo 



^1,0 



27r 



E 



1 
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^ + o + (M - ^O) 



2tt 



+ 



E 

1=0 



f3 



1 + :^] +^{p-£o) 



These series can be expressed in terms of the Hurwitz C- function |44l |45] . Taking 
into account that this function has branch cut discontinuities in the complex plane 
of its second argument running from to — oo, we get 
(6.11) 



E 

/ — — oo 



A 




-I- sign (pi - So) p 



2-K 



Im~^oI 



z-sign(/x-£o) P . _ ,\ 1 

The expression in braces has a vanishing analytic continuation at u = 0. As a 
consequence, the derivative of the right hand side at u = docs not depend on the 
arbitrary scale A. Then, the contribution of this Landau level to log DetD per unit 
area and at low temperature is given by |44[ 145] 
(6.12) 



oo 



du ^ 

^ — — oo 



■ 













1 



/3|M-fo| 



Taking into account that we are interested in the mean number of particles 
with respect to the neutral material (i.e., with all the Landau levels with negative 
energy filled, which corresponds to = 0), from Eq. (j6.2[) . we finally get for the 
contribution of the Landau level with energy to the Hall conductivity at zero 
temperature 
(6.13) 



B 

— a 

e 



xy 



Be d 



d 



= ^{ir [(^0 - m) e (£o - /i)] - ^ [(fo - Ai) e [Eo - /.)] 



27r I 9/i 



f^{-e(£o 

ZTT 



Be 
2^ 



{Q{EQ)Q{^i~EQ)-e\ 



Eo) 8 {Eo ~ fi) } 



Therefore, for positive Eo, we get a contribution (+|;;) to axy if /i > £o and zero 
otherwise. On the other hand, for negative Eq we get a contribution {~^) ii ^ < Eg 
and zero otherwise. 

At this point, it is clear that a similar result is obtained for any other Landau 
level. Indeed, the same considerations can be done for any other energy eigen- 
value E. Then, for a given chemical potential /z, there is only a finite number of 
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Figure 3. The Hall conductivity for a realistic value oi w = 10 



nonvanishing contributions and the Hall conductivity reduces to 



(6.14) 





where the sums run over the energy levels. Indeed, this is a rather general result 
for the Landau problem of a fermionic system. 

For the model of interest in the present paper, the spectrum to be considered 
is the union of those of and Hj^' , £ E Speci/x U SpeciJK' , which takes into 
account the states around both Fermi points. But they must be taken with an 
additional degeneracy corresponding to the two polarizations of the electron spin, 
which has played no role up to now in this analysis. 

In full units we get 



, 0<£<Ai 



2w 



'-). This is the characteristic of the 



anoma- 



(6.15) 

Notice that 
(6.16) 

for small (but > \£o\ = 

lous integer quantum Hall effect found in graphene [TH |T71 [THl HH] , which shows a 
nonvanishing Hall conductivity for small (positive or negative) Fermi level. 

In Figure 3, the Hall conductivity of this model as a function of ji (equal to the 
Fermi level at zero temperature) has been plotted for w = 10"^. Notice that it is 
not possible to appreciate the structure of each step as a double step (due to the 
small degeneracy breaking remarked in Eqs. (j5.2ip and (I5.22[) ) with such a realistic 
value of w (See Figure 4). Also notice that the Hall conductivity of our model 
vanishes in a small neighborhood of the origin, for |^| < 2£^IR_ 
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h 






cr 






1.6 




1.4 




1.2 




1.0 




0.8 




0.6 





1.405 1.410 1.415 1.420 1.425 VfVeS 

Figure 4. The structure of the step m the Hall conductivity from 
i to |, with w = 10^. 



7. Perturbation with the quadratic nearest neighbor interaction 

TERM 

Notice that the quadratic term in the nearest neighbor dispersion relation in Eq. 
(|5.8p can be obtained from the solutions of the free problem, Eqs. p. 21) and 
through the (rotational symmetry breaking) perturbation 

^ r. o \ — 1 



(7.1) 
Indeed, 

(7.2) 



AiJVk(x) = --«Fe*''-^fc2 [4fci'-k2] (k 



We will assume that, in the presence of an electromagnetic field A, derivatives 
change into covariant derivatives and the perturbation just turns into the Hcrmitian 
expression 
(7.3) 

a \ r 2 2I r 2I 

/SH = --VF \{p2 - eA2) 4(pi - eAi) - (p - eA) (p - eA) cr • (p - eA) 



+cr • (p - eA) (p - eA) 



4(pi - eAif - (p - eA)' (p2 - 6^2) 



Notice that there is an order indeterminacy in this definition. We will come back 
later on this point. 

We now consider a constant magnetic field as in Eq. (|4.ip . Then, 

(7.4) pi — eAi — pi — \/eB p , p2 — eA2 — P2 — eBxi — —VeB q , 
and AH reduces to 

(7.5) AH^Sj+Sj\ 
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where 
(7.6) 



o 



When S) is apphed to the functions in Eq. (|4.14p . taking into account that (See 
Eq. gS}) 



(7.7) 



[pai - qa2] ?/'fe,«,s(x) = < — Sn.s - 7; \P' 

w I eB 2 



where £"„ s is given in Eq. (|4.29l) , we get 
(7.8) 



■SoV'fc,«,s(x) = - WFveS 



W 



[2n + 3] 



2tt 

3(2n + 3)] Hn+i{q) 



\ 



1 - sy/l + 8w^{n + l) 



2w 



{f|£„,,-l(2n+l)} 
[2n + 1] 



[Aq^-3{2n + l)] H,M 



The first order perturbation on the eigenvalues of the Hamiltonian are given by 
the matrix elements 
(7.9) 

{lpk',n.s,^H'ipk,n,s) = {'(Pk',n,s,^'ipk.n,s) + i^'ipk',n.s,'ipk,n,s) 



dqe'^ q ■ 



^ 5{k ^ k') avpVeB 
{^£„^,-i(2n + 3)} 



[2n + 3] 



[4q^ - 3{2n + 3)] Hn+i\q) + 



+ 



Vl + 8w2(n + l) 
2^ 



{ag„,.-l(2n+l)} 
[2n + 1] 



[Aq^ -3{2n + l)]Hn\q) } ^ 



since it is the integral on the whole line of an odd function of q. 

As previously noticed, there is an order indeterminacy in the expression of Sj 
in Eq. (|7.6p . since the operators in the right hand side of Eq. ()7.ip commute. 
But any order one chooses will lead to an odd in q integrand in the expression of 
i'4'k' ,n,s, '^Hipk.n,s), giving then a vanishing result. 

It is also immediate to show that Sjipo = 0, with ipQ given in Eq. (|4.2ip . 

Therefore, we get no modification of the energy eigenvalues at first order in 
perturbation theory with AH in Eq. ()7.3p (or with any permutation of the factors 
in Sj in Eq. ()7.6p ). This means that, in the presence of a constant magnetic field 
orthogonal to the plane, the spectrum of the Hamiltonian of our model seems to be 
rather stable against the rotational symmetry breaking perturbation in Eq. (j7.5p . 

In particular, the Hall conductivity of the model described in Section [6] and the 
explanation of the anomalous integer quantum Hall effect based on it is not modified 
by the first order perturbation theory in AH, operator which incorporates the (next 
to leading) quadratic term in the nearest neighbor interactions. 
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8. Crossed Magnetic and electric fields 



Let us consider also an electric field in the direction of xi. So, we can take the 
electromagnetic field as 

(8.1) A -.= 3x162, Ao:=-Exi, 

also in the Coulomb gauge. 

In this case the Hamiltonian is 



H = 



(P - eA) 
2 m 



vp a ■ (p — eA) — e^o 



In the m — )■ 00 limit a pseudo-Lorentz transformation of the linear Hamiltonian 
allows to solve the eigenvalue equation for crossed electric and magnetic fields in 
terms of the solutions of the case in which there is only a magnetic field applied 
perpendicularly to the plane of graphene |46[ Wf\ H5]. But, in this nonrelativistic 
model, the quadratic mass term prevents to get the solution in this way. 

Translation invariance in the X2 direction suggests that the solutions are again 

Lp(xi) 



of the form 'I'(x) 



, and the eigenvalue equation reduces to 



Pi 



{eB) 



Xl 



k 

eB 



X + 2meExi > cp 



(8.3) 



'2mvF ipi+ leB ( ^ ^ ) r X , 



A + 2meExi > x 



-2mvF S Pi — ^eB xi 

eB 



with A = 2m£. 

It is convenient to define a new variable as q := ^feB — + 7w] ^^^^ 
correspondingly, p := — = af7! terms of which the eigenvalue equation 



writes as 
(8.4) 

where 
(8.5) 



+ Vil) = 2w <{ [p + iq] + 



2 I 2 

p +q 



E 
vpB 

E 



A := A 



mE 



mE 



- 2k 



with A = 2m£. 

For small electric fields, the last terms in the right hand sides of Eq. (|8.4p can 
be treated as the perturbation ^ ^^^ j f2 on a zero order Hamiltonian whose 
solutions are the same functions of q as in Section 21 with the equilibrium position 
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of the oscillator in the xi-direction displaced by an amount (— f^)- The zero- 
order energy eigenvalues are then those of Section 2] but shifted by an amount 
Elk - mE] 

On the other hand, since the perturbation is proportional to a2 , it has a vanishing 
mean value due to the orthogonality of Hermite functions (See Eq. (j4.14p ). Then, 
the spectrum is just rigidly shifted by the presence of an electric field perpendicular 
to the magnetic field, up to first order in perturbation theory in the small parameter 
(E/vfB). 

9. Conclusions and outlook 

In this article we have studied a simple nonrelativistic model, suggested by a non- 
Abelian magnetic quantization approach (based on a Heisenberg algebra deformed 
by the introduction of an external constant non- Abelian magnetic field) , to describe 
the low energy excitations of graphene. As previously mentioned, these states admit 
an effective description as excitations of a (2+l)-dimensional pseudo-relativistic 
field theory of massless Dirac fermions. 

Our effective model is based on a deformation of the Heisenberg algebra which 
makes the commutator of momenta proportional to the pseudo-spin. This can be 
interpreted as the commutator of two U{2) covariant derivatives in the case of a 
constant background non- Abelian magnetic field. The modification introduced in 
the Heisenberg algebra leads to wave functions of two components and a Hamilton- 
inan operator in the form of a 2 x 2 matrix, which reduces to the usually considered 
linear one in the large mass limit. 

We have solved the Landau problem for this planar system, with a constant and 
uniform magnetic field applied perpendicularly to the plane, explicitly obtaining 
the Hamiltonian (generalized) eigenfunctions and eigenvalues. We have pointed 
out that this construction can be directly extended to the case of negative values of 
the mass parameter. This is important since a negative mass allows to reproduce 
with our Hamiltonian the low energy expansion of the dispersion relation around 
the first Dirac point of graphene, K, when the leading order terms from both 
nearest and next-to-nearest neighbor interactions (linear and quadratic terms in 
the quasi-momentum expansion, respectively) are retained. 

The effective Hamiltonian for the excitations around the second Fermi point, 
Hk_' , is proposed to be obtained from Hiq through a transformation which maintains 
the Hamiltonian time-reversal times parity invariant. 

Through this transformation, the spectrum of the low energy states for the second 
Dirac point K' is simply obtained from that of K by a reflection around the origin. 
This leads to an almost doubly degenerate spectrum for quasi-particlcs and holes, 
where the degeneracy is broken by O {w^^) terms, with w — ~ 10'^ for 

realistic values of the parameters. Moreover, there is one quasi-particle state of 
energy £'q = eB/2\m\ and one hole state of energy £q = —eB/2\m\. 

The |m| — > oo limit reproduces the well known doubly degenerate spectrum of 
massless Dirac fermions, with one quasi-particle zero mode and one hole zero mode. 

The Hall conductivity of the model has been evaluated from the partition func- 
tion employing the ^-function approach to the associated functional determinant. 
This led to a rather general expression valid for the Landau problem of a fermionic 
theory. In the case under study, the result is consistent with the anomalous integer 
quantum Hall effect found in graphene, in which the (reduced) Hall conductivity as 
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a function of the Fermi energy takes half- integer values (See Figure 3). In partic- 
ular, Uxy — ±5 for small (but |ei?| > ^^^^ ) positive or negative Fermi level, 
respectively. 

Moreover, we have proposed an additional rotational symmetry breaking term, 
to be treated as a perturbation on the Hamiltonian, which reproduces the next- 
to-leading (quadratic) term in the nearest neighbor interaction contribution to the 
dispersion relation with no external field. With the perpendicular magnetic field 
present, we have shown that this term does not modify the spectrum at first order 
in perturbation theory. So, it neither changes the behavior of the Hall conductivity 
described above. 

Finally, we have shown that an electric field perpendicular to the magnetic field 
produces just a rigid shift of the spectrum at first order in perturbation theory. 

An interesting problem under study refers to the case in which an Aharonov- 
Bohm singular flux is also present, which can simulate a defect in the lattice. This 
will be the subject of a forthcoming paper. 
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Appendix A. Lagrangian, conserved current and generating 

FUNCTIONAL 

In this Appendix we construct the Lagrangian for our model and deduce the 
conserved current associated with its U{1) gauge symmetry. We also discuss the 
behavior of the generating functional in relation with the Hall conductivity. 

The system under consideration can be described by the Lagrangian 



C 



(A.l) 



-— I [(P -eA + Oct) i/-]^ • [(p - eA + Oa) ij;] - 2e'^tl)^\ 
2m I J 



The variation of the independent dynamical variables and 4' leads to the 
Euler-Lagrange equations in the usual way: 



(A.2) 



dC _ 



dt 



dC 



dC 



a(vv't) 



1 

2m 



(p - eA + eay - 20^ 



^ = 0, 



and similarly for ip. 



The Lagrangian in Eq. (jA.l[) has a U{1) gauge symmetry. Indeed, it remains 
invariant against the following local change in the dynamical variables and the 
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electromagnetic field: 

ip{x) e'-""'''''^ ij}{x) 5tp{x) = iea{x)ip{x) , 

(A.3) 7/;t(a;) -> V^^(x)e-*'="(^) ^ SiP^x) ^ ~iea{x)^l;\x) , 

Ai,{x) -> Ap(a;) + 9^a(x) . 

Then, the Noether's Theorem implies the existence of a locally conserved current 
given by 

Therefore, the charge density is 

(A.5) jO = eVV 

and the current density 

(A.6) j = |z (W^ ip-ip'' V^) + 2V/ (~eA + 6*0-) V| . 

2m 

It is a simple exercise to verify that j'^ is conserved as a consequence of the 
equations of motion, Eq. (|A.2p and its complex conjugate, 

(A.7) dtf-V-j = 0. 

The generating functional of (current) Green's functions of this two-dimensional 
system reads as 



, I dt d xCix) 

(A.8) Z[A]^e'^^^^:^ Vij^Vije 



where T[A\ is the effective action. In particular, we have 
,A.) 

and the second functional derivative gives minus the current-current correlation 
function. 

For a system with a gap in its spectrum (as is the case of our model with filled 
Landau levels), the effective action is a local gauge invariant functional of the 
electromagnetic external field [301 E] ^ 

(A.IO) T[A]^ j dt j d^xCeff{A^{x),d,A^{x),...) . 

In general, in an asymptotic weak field and weak gradient expansion, gauge in- 
variance determines that the effective Lagrangian Cej / appears as an expansion in 
terms of the field intensities E and B and its derivatives. But in 2-f 1 dimensions 
there is an additional possibility, since a Chern-Simons term can be the dominant 
one in this expansion (See, for example, [40l|41] and references therein), 

(A.ll) C,ff = ^ ^^""^A.d^A^ - \ E^p^^E, -lxB' + ..., 

where e^'^''^ is the completely anti-symmetric tensor with e°^^ = 1 and the ellipses 
represent higher order terms. 
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In fact, it is the Chern-Simons term which determines the HaU conductivity since 



The first equation in IA.I3I implies that the Hah conductivity a^y — 57 (or a^y — 
Ke'^/h in fuU units). Then, 



It can be shown that K takes integer values if the (many-body) ground state 
is not degenerate and has a finite energy gap |101 HU HH] ■ In this case, K has a 
topological character: It can be related to a Berry phase which, in turns, takes the 
value 2tt times an (integer) Chern number. 
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